Coexistence of excited states in confined Ising systems 



O 

o 
o 

(N 

< 
(N 



a 

e 

I 

a 
o 

(N 

On 

O 
O 

o 

a 

o 
o 



X 
J3 



Andrzej Drzewinski 
Institute of Low Temperature and Structure Research, Polish Academy of Sciences, 
P.O.Box 1410, 50-950 Wroclaw 2, Poland. 
(February 1, 2008) 

Using the density-matrix renormalization-group method we study the two-dimensional Ising 
model in strip geometry. This renormalization scheme enables us to consider the system up to 
the size 300 x 00 and study the influence of the bulk magnetic fleld on the system at full range 
of temperature. We have found out the crossover in the behavior of the correlation length on the 
line of coexistence of the excited states. A detailed study of scaling of this line is performed. Our 
numerical results support and specify previous conclusions by Abraham, Parry, and Upton based 
on the related bubble model. 

PACS numbers: 05.50.-fq, 68.35.Rh, 75.10.Hk 

The understanding of classical systems in confined ge- 
ometries has been a chahenge for several years [0-^. 
Among such investigated systems are fluids or magnets 
confined between parallel walls. Studies of finite-size ef- 
fects have not been limited only to the vicinity of the 
critical point, but also to the first-order phase transi- 
tions, which are less known. 

In this Raport we consider the two-dimensional Ising 
system on a square lattice in strip geometry {L is width 
of the strip) with the Hamiltonian 



A similar problem of higher-order corrections, but 
to the Kelvin equation (the scaling of the bulk coexis- 
tence field in the presence of the parallel surface fields) 
has been studied recently Using the density-matrix 
renormalization-group method (DMRG) §,0 it was 
found that for a large range of surface fields and temper- 
atures corrections are not compatible with the behavior 
predicted by the existing theory . It is one of reasons 
why we have checked out here the predictions given by 
the bubble model. 
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where the coupling J > 0, 77 is the bulk magnetic 
field and ai = ±1. The first sum runs over all nearest- 
neighbour pairs of sites while the second sum runs over 
all sites. 

Even such a simple model has an interesting crossover 
governed by the bulk magnetic field |^,|| , which value 
depends on temperature and the size of the system. The 
borderline Hx{T\L) divides two different L and H de- 
pendencies of the correlation length ^. Using the bubble 
model ||] Abraham et al. found ||] that at subcritical 
temperatures one has 

lli = P{T)L\Hl for < < (2) 

1/e = R{T) + ^(T)|i/|2/3, for \H\ > H^, (3) 

where P{T) = 2m/kBT, R{T) = 2ao/kBT. Here, 
S{T) is an unknown positive coefficient. Furthermore, 
m and ctq are the bulk spontaneous magnetization and 
the interfacial tension, respectively. 

The bubble model studies concluded that Hx{T;L) 
scales towards the first-order line according to the form 



(4) 



b) 



+m ( 


-m 


/ +m 




+m 


-m 





H^{L;T) « A{T)L°' + B{T)L'< + C{T)L^ 
where a = —1, 7 — —5/3, and 6 — —7/3. 



FIG. 1. The dominating configurations in strip geometry: 
a) Q<\H\< H^, b) \H\ > H^. 

Abraham et al. argued that the mentioned crossover 
occurs because the class of dominating configurations de- 
termining the behavior of correlation functions changes 
from a single connected loop for \H\ > Hj. to two dis- 
connected closed loops |i7| < (in cylinder geometry). 
In our case, where the free boundaries are present, for 
\H\ < Hx, the dominating configurations consist of suc- 
ceeding pieces of a strip with opposite magnetizations 
|]l2| . For \H\ > Hx the most important configurations 
contributing to the correlation function are again closed 
loops including domains of opposite magnetizations (see, 
Fig.l). 

In order to analyze this problem beyond the bubble 
model, we can use the transfer-matrix (TM) calculations 
p3[ . However, it is well known that to obtain satisfac- 
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tory finite-size scaling results, one should consider large 
enough systems JT^ . This may, in turn, complicate calcu- 
lations or even make them impossible. To overcome this 
problem we have applied the DMRG method for two- 
dimensional systems based on the TM approach. Pro- 
viding a very efficient algorithm for the construction of 
the effective transfer matrices for large L this method 
was successfully employed for a number of problems (for 
which no exact solutions are available, e.g. for nonvanish- 



ing bulk fields) |15-18|. Using it we were able to analyze 
the system in full range of temperatures and the bulk 
magnetic field for strips of widths up to L = 300. For a 
comprehensive review of background, achievements and 
limitations of DMRG, see Ref. 

We first calculated the free-energy levels 
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for z = 0, 1, 2, . . ., where Xi are the eigenvalues of the 
TM arranged in order of decreasing magnitude. Because 
the inverse (longitudinal) correlation length can be de- 
fined as 



l/^L) = log(Ao/Ai), 



(6) 



and the lowest free-energy level does not cross others, 
especially important are the values of the bulk magnetic 
field Hx{T;L), where the first and the second excited 
states cross each other. In such a case we can observe 
the crossover in the behavior of the correlation length. 

Let us first analyze the structure of the TM low-lying 
levels as a function of the bulk magnetic field H at fixed 
T. At very low temperature they should behave practi- 
cally in the same way as the ground state energy. There- 
fore, it is worthy first considering the ground state prop- 
erties of the system. 

Let us define the configuration of a row for the strip in 
the following way |cti, (72, cri_i, ctl), where the values 
of CTi are denoted ± for simplicity. At zero magnetic field 
H the two states with all spins positive | -I- -f ■ • ■ -I- or 

negative | • ■ • ) have the same energy. The extra 

magnetic field term splits both states and the energy per 
spin is 
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±H. 
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Assuming H > the | ++••• + +) state is always 
the singlet ground state. In order to find the first ex- 
cited states we have to fiip the first or the last column 
{i = 1, L) in the previous configurations. In this way we 

get the four states |— -I- • • ■-!-+), |++ • • ■ + — ), | + — • ) 

and I • • • — h). The magnetic field sphts this level 

into two doublets and for the two first states their energy 
decreases when the H increases according to the equation 



-J(2 - -) - H{1 



(8) 



Therefore, we expect the crossing of the singlet state 

) with the doublet |-H h-H), |-KH h-) 

at a certain value of the bulk magnetic field 



H4T = 0;L) = 



J 



(9) 



Note, that for T Eq.(|) reduces to Eq.(|) provided 
A{T) J and B{T),C{T) 0. 
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FIG. 2. The lines of coexistence of the excited states for 
different strip width L. Tc is the bulk critical point, whereas 
the thick solid line denotes the bulk first-order line. The ar- 
rows point at the inflection points where Hx {T; L) has a local 
minimum as a function of the temperature. T'{L) describes 
an end point where (T; L) ends on the H = axis. The 
dashed lines are as guides for eye. 

At finite temperatures we do not have any more real 
crossing points but, so-called "the regions of avoided level 
crossing". At = the first two levels are separated 
as /i — /o ^ exp(— (Toi/fcsr), so they are asymptoti- 
cally degenerated for _L — > cxd. The region of avoided 
level crossing continues for nonzero magnetic fields up to 
\H\ ~ exp{-~aoL/kBT). When we are interested in the 
behavior of ^, we have to consider the second and third 
eigenvalues of the transfer matrix , where asymptotic 
degeneracy is also present for /i and /2 ||2^ . It is assumed 
that the difference /2 — /i and the avoided level cross- 
ing region centered on are of order exp(— CL), where 
the coefficient C may be H and T dependent. It makes 
sense to discuss algebraic shift of the value Hx{T; L) for 
L ^ oo. In order to find out the value of in finite tem- 
peratures at fixed L, we identify with the value of the 
bulk magnetic field where the second free-energy level /i 
has a maximum and the separation /2 — fi is the smalest 
one. The curves Hx{T\ L) present the coexistence of the 
excited states and are shown in Fig. 2. 
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The curves indicate the phase boundaries between the 
two phases with different dependencies of ^ on L and H . 
As i — > cxo the coexistence hues of the excited states 
shift towards the H = axis that is intuitively clear at 
T = 0. Since the width of the strip increases the energy of 
configurations where only one column of spins is flipped 
(Eq.(||)) are close and close to the energy of configura- 
tions with all spins pointed in one direction (Eq.(|^)), so 
in the i — > oo limit one has Hx = 0. 

Let us discuss the scaling of the coexistence line 
Hx{T;L) to the bulk first-order line (Fig.(|)). To ver- 
ify the bubble model predictions (Eq.(||)) we have cal- 
culated series of values of Hx{T = const; L) for L — 
20, 40, . . . , 200 and for temperatures ranging from T « 
0.44Tc up to T w 0.99Tc. 

Table. Scaling exponents of Hx{T;L) 
extrapolated from the DMRG data. 
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-0.67(2) 


2.20 


-1.002(6) 


-0.69(1) 


-0.68(3) 



Table shows the values of scaling exponents obtained 
from the DMRG data. Using the powerful extrapolation 
technique, the Bulirsch and Stoer (BST) method |2^, we 
obtained an excellent agreement with Abraham et ai. 

In order to get the A coefficient in Eq. (^) one can com- 
pare Eqs. (||) and (^). They have to agree at the value 
H — in the termodynamic limit, which implies the 
following relation M: 



A{T)^aoiT)/m{T). 



(10) 



In Fig. II our data reconstruct this curve very well. To 
the best of our knowledge, the coefficients B{T) and C{T) 
in Eq.(|^) have not been yet determinated, but our numer- 
ical results can predict their temperature behavior. 

Close to Tc the validity of Eq. (^) is limited because the 
scaling of points of the Hx curve is governed by the bulk 
critical point. In order to study it in detail we have con- 
sidered characteristic points of the upper part of the Hx 
curve: the inflection points {Hc{L) ,Tc{L)) and the end 
points T'{L) (see, Fig.(|^)), where the following scaling is 
expected: 



T,{L) = {T,-T,{L))/T, 
H,{L)^L-y". 



(11) 



Here, i/t — I and yn ~ 15/8 are the thermal and mag- 
netic exponents of the two-dimensional Ising model. 

To verify the scaling to the critical point {H = 0,r = 0) 
we found out the inflection points and the end points for 



L = 30, 60, 100, 130, 160 and 200 using subsequently the 
BST technique. We have examined the scaling form ( |Tl| ) 
for L — > oo and found very good agreement 

Tc ^ 0.00006(6), and yr = 1.005(5) 
He = -0.0006(6), and yn = 1.876(8) 
t' = 0.0000(3), and yr = 1.006(6) 

Note, that for T and L ^ oo we can repro- 

duce the scaling form ([l|) from Eq. (|^) by assuming that 
A{T) 0, B{T) 0, and C{T) oo. This is in 
agreement with our numerical estimations for scaling co- 
efficients as depicted in Fig.^. Of course, this relation is 
not valid at Tr. 
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FIG. 3. The coefficients of the scaling of the coexistence 
line Hx{L;T) to the bulk first-order line. Solid line denotes 
the analytical result determined in the bubble model. The 
symbols describe our numerical results. The dashed lines are 
as guides for eye. 

In order to analyze the behavior of the correlation 
length we have derived 1/^ for L between 100 and 300 
in temperatures below Tc{L). To examine the form of 
Eq. (^ firstly we have confirmed the linear dependence of 
the coefficient on L. Next we compared our numerical 
resuhs with the coefficients P{T) and R(T) in Eqs.(|,|l). 
What is more, we presented the temperature dependence 
of the S{T) coefficient which was not determined in the 
bubble model (see, Fig.^). 

When temperature raises, more and more complex con- 
figurations on the Ising strip contribute to the free en- 
ergy in contrast to the assumption of the bubble model 
Consequently, in high temperatures the validity of 
Eqs. (^,||) is limited to a narrow range of H. The bub- 
ble model predictions are also spoiled by the presence of 
strong bulk magnetic field. That is why, the higher is 
temperature the smaller H is necessary to recover the 
linear dependence of 1/^ on H, as in Eq.(|). Similarly, 
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when T Tc{L) the regime with the i?^/^ dependence 
of 1/^ (Eq.(^), close to the right side of the coexistence 
hne, shrinks to zero. 



critical reading of the manuscript. This work was sup- 
ported by the Polish Science Committee (KBN) under 
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FIG. 4. The coefficients of the correlation length in 
Eqs.(^,^. Solid lines denote the bubble model results: 
P(T) = 2m/kBT and R{T) = 2ao/kBT. The symbols de- 
scribe the corresponding DMRG results. 

In conclusion, we have used the density-matrix 
renomalization-group method to obtain reliable informa- 
tion about the two-dimensional Ising model in the bulk 
magnetic field. We have confirmed the crossover related 
to the correlation length analyzed before for the bubble 
model 1^ . Our study has not been limited to subcritical 
temperatures and small bulk fields. We have confirmed 
Abraham et al. predictions for the scaling of the first- 
order line in the subcritical region. Morover, we have 
established the precise scaling form for the bulk mag- 
netic field by numerically determining coefficients B{T) 
and C{T) in Eq.(^). Furthermore, we have extended the 
analysis of the bubble model to critical region verifying 
that the scaling behavior is governed by the bulk critical 
point. Finally, we numerically confirmed the magnetic 
field dependence of the correlation length, simultaneosly 
extracting previously unknown coefficient S(T) in Eq.(^. 
Above results demonstrate that for two-dimensional clas- 
sical systems the DMRG technique provides significantly 
accurate data for studying equilibrium properties of large 
systems in a nonvanishing bulk magnetic field. 
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